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Abstract: 

The new result for the third-order QCD corrections to Rf,+^-, unhke the old, incorrect result, 
is nicely compatible with the principle-of-minimal-sensitivity optimization method. Moreover, 
it leads to infrared fixed-point behaviour: the optimized couplant, Os/vr, for i?e+e- does not 
diverge at low energies, but "freezes" to a value 0.26 below about 300 MeV. This provides some 
direct theoretical evidence, purely from perturbation theory, for the "freezing" of the couplant 
- an idea that has long been a popular and successful phenomenological hypothesis. We use 
the "smearing" method of Poggio, Quinn, and Weinberg to compare the resulting theoretical 
prediction for R^+e- with experimental data down to the lowest energies, and find excellent 
agreement. 



1 Introduction 



The calculation [|^, ^ of the third-order (next-to-next-to-leading order) QCD corrections to 

Re+e- ■ 

i?g+g- = (Tfot(e^e~ — > hadrons)/(T(e'''e~ fi~^fi~), (1-1) 

provides valuable empirical information on the behaviour of perturbation theory in QCD. This 
paper is concerned with "optimized perturbation theory" (OPT) and is motivated by three 
questions which the Re+e- calculation can answer: |Q] 

(1) Does perturbation theory seem to be well behaved? Is the third-order "optimized" result 
in reasonable agreement with the second-order "optimized" result? What can we learn about 
the error estimate? 

(2) Is the optimized couplant, a = as /it, smaller in third order than in second? The "induced- 
convergence" picture suggests that the optimized couplant a^"^ as determined by the nth- 
order optimization equations, will tend to decrease as the order n increases. In this way "op- 
timization" could lead to a convergent sequence of perturbative approximations, even if the 
perturbation series in any fixed renormalization scheme is divergent [|| . 

(3) Does one find infrared fixed-point behaviour? A third-order calculation is a prerequisite for 
addressing this question in "optimized" perturbation theory, and the answer basically depends 
upon whether the invariant p2 (defined below) is negative or not [^. 

It is striking that, with the originally published third-order result [0, the answer to all three 
questions was "No" — while, with the new, corrected, result the answer to all three questions 
is "Yes" . The purpose of this paper is to elaborate on these three points, and especially to discuss 
the infrared fixed-point behaviour |^. We do not share the pessimistic attitude of Chyla et al 
10, ^ to the infrared results. If one believes in OPT, the infrared results - though quantitatively 
uncertain - are qualitatively unequivocal: We propose to take them at face value and compare 
them to experimental data ||8|. 

The plan of this paper is as follows: Sect. 2 reviews OPT, and applies it to Re+e- third- 
order, with particular emphasis on the infrared limit. Sect. 3 compares the predicted Re+e- 
with experimental data, using Poggio-Quinn- Weinberg (PQW) smearing |T^. Sect. 4 briefly 
discusses the phenomenology of a "frozen" couplant. Conclusions are summarized in Sect. 5. 
Some technical matters are relgated to the appendices. 

2 Optimized Pertubation Theory and Fixed-Point Behaviour 
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2.1 The principle of minimal sensitivity 

We begin with a few words about the principle of minimal sensitivity, upon which OPT is 
based. It deals with any situation where an exact result is known to be independent of certain 
variables, but where the corresponding approximate result depends upon those variables, and 
hence is ambiguous. (In the QCD context, physical quantities are Renormalization Group 
(RG) invariant but perturbative approximations to them are not, due to truncation of the 
perturbation series.) The philosophy is that such a non-invariant approximant is most believable 
where it is least sensitive to small variations in the extraneous variables, because this is where 
it best approximates the exact result's vital property of being completely insensitive to the 
extraneous variables. 

A simple example is perhaps the best way to convey this idea. Consider the quantum- 
mechanical problem of computing the eigenvalues, £'fc, of the quartic-oscillator Hamiltonian: 

+ ^3.4^ (2.1) 



where [x,p\ = i. Suppose we do standard perturbation theory, but with |T^: 

Ho = i(/ + O^x^), Hint = Ax^ - iJ^^x^. (2.2) 

This introduces an "extraneous variable" fi, and the approximate eigenvalues so calculated will 
be f]-dependent. For example, first-order for the kth. eigenvalue gives: 



4'^ = + m + |^(2fe' + 2k + 1). (2.3) 



However, we know that the exact eigenvalues are O-independent. Therefore, it is sensible to 
choose so that the approximant, Ej^\ is minimally sensitive to i.e., 



■ {2k^ + 2k + l)' ' 
oA 



+ 



(2.4) 



(Quite generally, we shall use an overbar to denote an "optimized" value.) This gives the 
"optimized" result: 



4^^(opt.) = -(A; + i) 



■ (2A;2 + 2A; + 1 ' ' 
oA 



(k + h) 



(2.5) 



This simple formula fits the ground-state energy to 2%, and all other energy levels to within 
1%. The secret of this success is the "optimal" choice of il, which is different for different levels. 

One may proceed to the calculation of higher-order corrections for some specific eigenvalue 
{e.g., the ground state, k = 0). For any fixed 0, the perturbation series would diverge, but if 0, 
is chosen in each order according to the "minimal sensitivity" criterion (which gives an that 
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gradually increases with order), one finds quite nice convergence [12j. This is an example of the 
"induced convergence" mechanism 13]. 

One may also use the same method to obtain accurate approximate wavef unctions, ipkix), 



from first-order perturbation theory [14|. Here the "optimal" will be a function of x; in 
particular, it will be proportional to |x| at large \x\, thereby converting the Gaussian dependence 
exp(— ijlx^) into the correct large-|x| behaviour. A variety of other examples and applications 
can be found in Refs. HII, |l| . Some examples of QCD applications can be found in Refs. 



16, 



2.2 RG invariance and optimization 

We next review "optimized perturbation theory" (OPT) |^ as applied to the QCD corrections 
to the R(,+e- ratio. Ignoring quark masses for the present, we may write R^+e- = Q^iX + T^)^ 
where TZ has the form: 

7^ = a(l + rla + r2a^ + ...), (2.6) 

and depends upon a single kinematic variable, Q, the cm energy. OPT is based on the funda- 
mental notion of RG invariance , which means that a physical quantity is independent of the 
renormalization scheme (RS). Symbolically, we can express this by: 

(i7^ _ 57^ da dTZ 

~ d{RS) ~ d{RS)^ d(RS)~d^' ^'^ 

where the total derivative is separated into two pieces corresponding to RS dependence from the 
series coefficients, rj, and from the couplant, respectively. A particular case of Eq. is the 

familiar equation expressing the renormalization-scale independence of TZ: 

''|+''<'"|;)'^ = »' (^^«> 

where 

f3(a) = = -ba^n +ca + 020^ + ...). (2.9) 

The first two coefficients of the /3 function are RS invariant and, in QCD with A'^^ massless 
flavours, are given by: 

(33 - 2Nf) _ 153 - 19A^/ 



6 2(33 - 2Nf) 

When integrated, the /3-function equation can be written as: 



(2.10) 



da' ^ fi' dii' , , , 

+ C= I -^i = ln(^/A). (2.11) 



P{a') J A /i 
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where C is a suitably infinite constant and A is a constant with dimensions of mass. The 
particular definition of A that we use corresponds to choosing 

da' 



C 



6a'2(l + ca') 



(2.12) 



(where it to be understood that the integrands on the left of (2.11) are to be combined before 
the bottom limit is taken). This A parameter is related to the traditional definition [18| by an 
RS-invariant, but A'^y-dependent factor: 

ln(A/A) = (c/6) ln(2c/6). (2.13) 

The A parameter is scheme dependent, but the A's of different schemes can be related exactly 
by a 1-loop calculation 20|. As is usual, we shall regard A^ (for 4 flavours) as the free 
parameter of QCD [21|. 

From Eq. (2.11) it is clear that a depends on RS only through the variables /i/A and 
C2,C3,..., the scheme-dependent /?- function coefficients. The coefficients of IZ can depend on 
RS only through these same variables, because of RG invariance, Eq. ( |2.7D . Therefore, these 
variables provide a complete RS parametrization, as far as physical quantities are concerned [^. 
Thus, we may write: 

a = a(RS) = a(r, C2, C3, . . .), (2-14) 

where 

T = hhl{^l/K). (2.15) 

The T variable is convenient and also serves to emphasize the very important point that RS 
dependence involves only the ratio of /i to A. "Optimization" does not determine an "optimum 
/i", but it will determine an optimum r. 

The dependence of a on the set of RS parameters r and Cj is most easily obtained [22| by 
taking partial derivatives of Eq. ( |2.11| ), varying one parameter while holding the others constant. 
This yields: 

(2.16) 



da 
dcj 



Pj{a) = -b(5{a) / dx 



2 ■ 



(2.17) 



Note that the I3j functions begin at order a^^^. 



The symbolic RG-invariance equation (2.7) can now be written out explicitly as the following 
set of equations: 



d_ 
d_ 



+ 



f3{a) a 



+ Pj{a) 



da 
d_ 
da 



7^ = 0, 

\n = o 



0=2,3,...). 



(2.18) 
(2.19) 
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These equations determine how the coefficients of TZ must depend on the RS variables. Thus, 
ri depends on r only, while r2 depends on r and C2 only, etc., with 

|1 = 1, (2,20) 

dro dro 

etc.. Upon integration one will obtain rj = /(r, C2, . . . , Cj) + const., where / is a known func- 
tion and the constant of integration is an RS invariant. Thus, certain combinations of series 
coefficients and RS parameters: 

p^(Q) = r-ri, (2.22) 

P2 = r-2 + C2 - (ri + ic)^ (2.23) 

etc., are RS invariant |3|, In the e~^e~ case, pi is a function of the cm energy Q, while p2 
and the higher-order invariants are pure numbers, dependent only on the number of flavours, 
Nf. 

Although the exact TZ is RG-invariant, the truncation of the perturbation series spoils this 
invariance. The nth-order approximant T?.*-""^ , defined by truncating TZ and the /? function to only 
n terms, depends on the RS variables r, . . . , c„_i. OPT corresponds to choosing an "optimal" 
RS in which the approximant TZ^"'^ is stationary with respect to RS variations; i.e., the RS in 
which 7^(") exactly satisfies the RG-invariance equations, ( 2.18 , 2.19 ). (Note that only the first 
(n — 1) equations will be nontrivial in nth order.) 

The second order approximant is 

TZ'^^^ = a{l + ria), (2.24) 

where a here is short for a^'^\ the solution to (|2.11| ) with /? truncated at second order. TZ^"^^ 
depends on RS only through the variable r. The optimization equation, from ( p. 18 ), is 

a'^ -a'^{l + ca)il + 2ria) =0. (2.25) 

This equation, together with the pi definition and the second-order integrated /^-function equa- 
tion, ( p. 11 ), uniquely determines the optimized result. (For details, see P, |l^].) 
The third order approximant is: 

=a(l + ria-hr2a2), (2.26) 

where now a is short for a^^\ the solution to ( |2.11 ) with /? truncated at third order. TZ^^^ depends 
on RS through two parameters r and C2, so there are two optimization equations, coming from 
(|2ll , [2l9|) . These can be reduced to g: 

(3f2 + 2fic + C2) + (3f2C + 2fiC2)a + 3f2C2a2 = 0, (2.27) 
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1(1 + (c + 2fi)a) - a = 0, where I = [ — 

Jo (1 

This integral can be done analytically, and is given by: 



dx 



(1 + CX + 02X^)2 



a(c^ — 2c2 + cc2a) 



(1 + ca + C2a^ 



4c2/(a,C2) 



with 



/(a,C2) 



2A 



In 



l + ia(c + A) 
l + ia(c-A) 



(2.28) 



(2.29) 



(2.30) 



where = — 4c2. (This assumes that A^ > 0, which proves to be true here.) The procedure 
for solving these optimization equations is discussed further in subsection 2.4, but we next 
discuss the infrared limit. 



2.3 Infrared limit and fixed-point behaviour 

Suppose we consider TZ at lower and lower cm energy, Q. Since c is positive for Nf < 8, 
the second-order (3 function has no non-trivial zero. Thus, in any RS, the couplant a*-^^ and 
approximant TZ^"^^ must become singular at some Q of order Aj^. In third order this may or 
may not happen, depending on whether the RS has a positive or negative C2. If C2 is negative 
then the couplant remains finite and tends to a "fixed-point" value, a*, which is the non-trivial 
zero of the third-order /3 function; i.e., the positive root of 



1 + ca* + C2a 



*2 



0. 



(2.31) 



Since fixed-point behaviour hinges on C2, which is scheme dependent, it is vital to have a sensible 
choice of RS In OPT the optimal C2 is determined by the optimization equations, and 
depends on Q somewhat. If C2 is negative as Q — > 0, then OPT will give fixed-point behaviour. 
The infrared limit of the optimization process was analyzed in Ref. and we briefly review the 
relevant results. 

Since /? vanishes at a fixed point, the r optimization equation, corresponding to ( |2.18| ), 
reduces to 

1 + (2fi + c)a* = 0. (2.32) 
Then, just by differentiating ( p.31| ) with respect to C2, one obtains: 



lim^4^^(a) 



da*_ 
dc2 



-,*2 



(2.33) 



(c + 2c2a*) 

(This can also be obtained, more laboriously, as the limit of ( p.l7| ).) Thus, the C2 optimization 
equation, corresponding to ( 2.19| ), becomes 



a + 



;i + 2fia* + 3f2a 
(c + 2c2a*) 



*2\ 



0. 



(2.34) 
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The two optimization equations yield: 



1 (1 + ca*) _ 2_ 
n = --^^, r, = --c,. (2.35) 



Using the expression for the invariant p2-, (p.23|), one obtains: 



C2 = 3(p2 + (2.36) 



Finally, substituting into the fixed-point condition ( 2.31| ), one finds Q: 

\ + ca* + 3p2a*^ = 0, (2.37) 

which determines a* in terms of the RS-invariant quantities c and p2- A positive a* exists if p2 
is negative, and the more negative p2 is, the smaller a* will be. 

2.4 Implementing the optimization procedure 

Returning to finite Q, we now consider how to obtain the third-order optimized approximant 
j^i^) numerically as a function of Q. As input, we need the values of pi and p2- Being invariants, 
they can be obtained from calculations performed in any computationally convenient RS. The 
calculations in the literature have used the "modified minimal subtraction" (MS) convention, 
with the renormalization point p chosen to be Q. The TZ coefficients are [p^ , |^: 

ri(MS; m=Q) = 1.9857 - 0.1153iV/, (2.38) 

r2(MS; m=q) = -6.6368 - 1.20017V/ - 0.0052iV| - 1.2395 (E%)7(3E9r^)- (2-39) 
The RS parameters of the MS{p = Q) scheme are: 

r(MS; /.=q) = 61n(Q/Aj^) = 61n(Q/Aj^) + cln(2c/6), (2.40) 



C2(ms) 



2857 5033 325 ,,2 

Nf H Nj 

2 18 54 



(2.41) 



16 (33 - 27V/) 

(The latter was first calculated in Ref. [^5|, and has recently been confirmed independently p^.) 
Substitution of these results into ( 2.22 , 2.23| ) gives the invariants. One can see explicitly that 



pi depends logarithmically on the cm energy Q, and on the free parameter of QCD, Aj^ 
However, p2 depends only on Nf. Since p2 turns out to be negative, one will find fixed-point 



behaviour in the "optimum" scheme [27, 28 1. Table 1 gives the fixed-point couplant values, 
determined from Eq. ( |2.37| ), for various Nj values [29|. 

Consider a world with Nf massless quarks, ignoring complications due to quark thresholds for 
the present. For simplicity we assume that the value of Aj^ is given, and our numerical results 
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use = 230 MeV for four flavours ||30[. For any chosen Q we then have definite numerical 
values for the invariant quantities pi, p2 (and 6, c). We need to solve for the optimum couplant, 
a, and the optimized coefficents, n, f2, and this will involve determining the RS parameters f, 
C2 of the optimal RS. These five variables are related by five equations; the two optimization 
equations, ( 2.27 , p.28 ), the pi, p2 equations, ( 2.22 , 2.23 ), and the integrated /3-function equation, 



(2.11), which for the /3-function truncated at third order becomes explicitly: 

r = - + cln(ca) - ^cln(l + ca + c^c?) - {(? - 2c2)/(a, C2), (2.42) 
a 2 



where f{a,C2) is given by ( 2.30| ). By substituting this last equation into the pi equation we can 



obtain fi explicitly as a function of a, 02- We can then rearrange the p2 equation to give f2 
explicitly as a function of a, C2- This leaves a, C2 to be solved for from the two optimization 
equations. Starting from an initial guess for a and C2, our procedure was to solve ( p. 28 ) numer- 



ically for a new a; and then (with the new a) to solve ( |2.27 ) for a new C2. We then iterated 



this procedure until the difference between successive solutions reached a specified tolerance. 
Further details are given in Appendix A. 

At very low Q we encountered technical problems with slow convergence of the iteration 
scheme. These are discussed in Appendix A. Nevertheless, with care it was possible to obtain 
accurate solutions at low energies. In Fig. 1 we show the optimized solution in the a, C2 plane as it 
smoothly approaches the fixed-point solution, which lies on the infrared boundary l+ca+C2a? = 
0. The figure shows two cases, Nf = 3 and Nf = 2. (In the real-world case we must switch 
from 3 flavours to 2 when we cross the strange-quark threshold. This requires a matching of A 
parameters, as discussed in Appendix B.) 

The optimized couplant a is shown as a function of Q in Fig. 2. Note that the effective 
couplant below 300 MeV is nearly constant at about 0.263, which is the Nf = 2 fixed-point 
value. Fig. 2 also shows the second- and third-order optimized results for TZ. The second-order 
result diverges at Q ~ 400 MeV, where pi{Q) vanishes. However, TZ^^^ remains finite, rising 
only to 0.33 at Q = 0. 

2.5 Illustrative results 

We pause for a moment to consider a comparison between the second- and third-order optimized 
results at moderately high Q. This exercise was performed by several authors fs^, ^ when the 
'old' third-order result Q was first published, and the results were disquieting. However, the 
new result [0] has transformed the situation, which is now very satisfactory. In Table 2 we give 



details for the two illustrative cases considered by Maxwell and Nicholls |31], namely Nf = 5, 
Q = 34 GeV, with either A,^ = 100 MeV or A,^ = 500 MeV. [Note, though, that the results 



8 



depend only on the ratio of Q to A^.] From Table 2 we see that between second and third 
order the optimized prediction TZ decreases only a few percent. With the 'old' result there 



had appeared to be a disconcertingly large increase |31, The new situation is much more 
satisfactory in other ways, too: In both examples the coefficient f2 now has a more reasonable 
magnitude, and the f\ coefficient has not changed so drastically from second to third order. 

The optimized couplant a now shows a marked decrease from second to third order. This 
is just what one would expect in the "induced convergence" picture of Ref. Q. In that pic- 
ture "optimization" induces convergence through a mechanism in which the effective expansion 
parameter, a, shrinks from one order to the next. Note that the 'old' result gave the opposite 
behaviour, with a apparently increasing from second to third order. 

The results also shed some light on the error-estimation question: If we knew just ^^^^ = 
a(l + fia), how might we estimate the error? Two estimates suggest themselves: (i) na^, where 
n is an order-one number, which presumes a well-behaved converging series, or (ii) |fia^|, the 
magnitude of the last calculated term, which is a typical error estimate for an asymptotic series. 
Knowing VS'^^ we can, presumably, get a much better idea of the actual error in 

7^,(2) from the 

difference 5 = 'R.^'^^ — TZ^"^^ . We have compared 6 with estimates (i) and (ii) over a wide range of 
Q/Ay^ values. 

Estimate (i), if we had assumed n ~ 1 or 2, would have been rather too optimistic. In 
fact, \6\ is between 7 and 14 times (for Nf = 4 and Q/Aj^ > 5). This is directly related 
to the size of the invariant p2 (which is about —14 for 4 flavours). [One can show analytically 
that 6 = p2a? + 0{a^) in the large-Q limit.] Of course, we could not know p2 until a third- 
order calculation was done. Arguably, though, 14 can still be considered an "order-one" number, 
especially in a theory that naturally involves numbers such as 4 (flavours), 3 (colours), 8 (gluons), 
etc.. 

Estimate (ii), based on the last calculated term, agrees with \5\ to within a factor of two 
either way for Q/Ay^ between 5 and 1000. At higher Q/Ay^ values this estimate would be 
overly pessimistic. However, we think that for present energies the estimate (ii) is perhaps the 
safest way to estimate the error. We suggest that it be used in QCD applications where only 
second-order results are known. 



2.6 Credibility of the infrared results 

We have stressed that OPT yields finite results for TZ down to Q = 0. The crucial question is, 
of course: How meaningful are these results? We would like to explain why, in contrast to other 
authors 0, ^, we take a positive attitude on this issue. 

Firstly, suppose we adopt the philosophy that the last calculated term in the "optimized" 
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perturbation series is a measure of the error. As we saw in this last section, this proved to be 
reasonable in the second-order case. In third order this gives \f2a^\ as the error estimate. Since 
TZ a, this implies a fractional error of 1^2^^^ I- In Table 3 we show some illustrative results 
at low energies, together with their estimated error. From this one can see that the behaviour 
of the series is quite satisfactory above Q = 1 GeV. The situation undoubtedly deteriorates at 
lower energies; by the time we reach Q = we have a series of the form 0.26(1 — 0.76 + 1.01) in 
which the higher-order terms are comparable to the leading term. While this is hardly a good 
situation, it is not completely disastrous; the corrections alternate in sign, and they do not dwarf 
the leading term. We believe that our error estimate, which grows to 100% at Q = 0, is not 
unreasonable: the result may well be off by a factor of 2, but is unlikely to be off by an order of 
magnitude. The qualitative conclusion, that IZ remains small (say, 0.3 it 0.3) at low energies is 
hard to escape. 

Secondly, it is instructive to view the use of QCD perturbation theory in the infrared limit as 
an extrapolation away from A'^^ = 33/2 [^4|. At Nj = 33/2 the leading /3-function coefficient, 5, 
vanishes (and hence c goes to — oo). For Nf = 33/2 — e, with e small and positive, there must be 
an infrared fixed point at a* ~ — 1/c = 0{e) |3^]. Perturbative calculations, even in the infrared, 
should then be meaningful if e is sufficiently small. Furthermore, one could naturally expect that 
the more orders in perturbation theory one has, the further one can extrapolate from Nf = 33/2. 
With sufficient orders one should be able to get infrared results down to Nf = 0, unless there is 
some unknown reason for the behaviour of the theory to change fundamentally at some critical 
Nf between 33/2 and 0. What does happen? Well, at second order, of course, one finds fixed- 
point behaviour, with a* = —1/c, provided c is negative, which requires Nf > 153/19 ~ 8, 
though Nf needs to be still larger if a* is to be reasonably small. In third order our results 
imply that, in the i?e+e- case, fixed-point behaviour — with moderately small a* values — does 
extend to Nf = 0. 

In the e ^ limit, a* tends to —1/c, and hence to (8/321)e. The small coefficient sug- 
gests that the natural expansion parameter of an extrapolation from Nf = 33/2 is not e but 
approximately e/40. One can verify that the third-order OPT results smoothly approach the 
limiting form as e ^ 0. For Nf = 16 one has p2 = —1724.4, and one gets a series of the form 
0.012(1 — 0.03 -|- 0.04). As Nf decreases, the behaviour of the series deteriorates, but it does so 
quite steadily; there is no dramatic change around Nf = 8 or any other Nf. 

In conclusion, our view is that the Nf = 2 infrared results, while quantitatively uncertain, 
are qualitatively credible. Having made this case in theoretical terms, let us now see what 
experiment has to say. 
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3 Comparing Theory to Experiment 



3.1 Re+e- including quark masses 

In this section we construct the theoretical prediction for R^+e- (allowing for quark masses) and 
discuss its comparison with experiment using the PQW smearing method. We limit ourselves 
to the region below 6 GeV, and we shall be particularly interested in the region below 1 GeV. 
To allow for quark masses in R^+e-^ we used the following approximate formula p 



Re+e- =3T.'l'^ Tiv,)[l + givi)n], (3.43) 

i 

where the sum is over all quark flavours that are above threshold (i.e., whose masses, rrii, are 
less than Q/2), and 

V, = {l-Am^/Q^)-2, 
T{v) = v{3-v'^)/2, (3.44) 



vr (3 + f ) /vr 3 



_2v 4 V2 47r, 

The coefficient T(f j) is the parton-model mass dependence and g{vi) is a convenient approximate 



form for the mass dependence of the leading-order QCD correction |10, 35 1. The higher-order 
corrections have been calculated only for massless quarks, so we simply evaluate TZ with Nf 
equal to the the number of above-threshold flavours. 

In our numerical results we used standard values for the current-quark masses = 5.6 

MeV, rud = 9.9 MeV, = 199 MeV, rric = 1.35 GeV. For Aj^ we used a 4-flavour value of 
230 MeV above charm threshold {Q > 2mc). Then, each time a flavour threshold was crossed 
as we decreased Q, we reduced Nf by 1 and computed the new parameter appropriate to 
the new Nf. The matching of A's is discussed in Appendix B. 

In this way we obtained the "raw" theoretical prediction for R^+e- shown in Fig. 3. For 
comparison, the flgure also shows the parton-model result {i.e., with the QCD correction term 
TZ set to zero). 

3.2 PQW smearing 

A direct comparison of the theoretical prediction with the experimental data is not possible, 
because there is no direct correspondence between the perturbative quark-antiquark thresholds 
and the hadronic thresholds and resonances of the data. However, a meaningful comparison is 



possible if some kind of "smearing" procedure is used |10, 37 1. We used the smearing method 
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of Poggio, Quinn, and Weinberg (PQW) who define the "smeared" quantity: 

R,^M^) = -j^ ds (,._Q2)2 + A2 - (3-45) 



In terms of the vacuum-polarization amphtude H, one can write Rpqw as ||T( 

2iRpQ^{Q; A) = n(Q2 + ^A) - U{Q'^ - iA). (3.46) 

In the hmit A — > this reduces to 2ii?g+e- j which is the discontinuity of 11 across its cut. 
However, a finite A keeps one away from the infrared singularities and nonperturbative effects 
that lurk close to the cut. The idea is to apply this smearing to both the theoretical and 
experimental i?e+e-'s and then compare them. 

In principle, the more orders in perturbation theory one has, the smaller one can take A 
However, this requires the full mass dependence of the higher-order corrections, which we do 
not know. In their leading-order study of the charm-threshold region, PQW used a value A = 3 
GeV^, and we shall use values of the same order of magnitude. We take a pragmatic view: the 
best choice of A is the smallest value that will smooth out any rapid variations in either the 
experimental or the theoretical R^+f,-. It turns out that this depends upon the energy region 
one is interested in. Around charm threshold a A of 3 GeV^ or more is necessary, while in the 
lowest energy region a A as small as 1 GeV'^ can be used. 



The integral in Eq. ( 3.45 ) was evaluated by numerical integration, after first making a 
change of variables s' — = Atan0. The computer routine was designed to take an input 
R-e+e- ) specified over a range to Qmax and to evaluate the integral over this range. A term 
was then added to account for the contribution from Qmax to oo, assuming that R^+e- remained 
constant above Qmax- The accuracy of the numerical-integration routine was tested against 
analytic results for several simple input functions. 

3.3 Experimental data and resonances 

The experimental data we used comes from a variety of sources: e+e~ vr+vr^ data in the p 
region and above from the OLYA/CMD and DM2 collaborations pi, |9|; Adone 772 data from 



1.4 to 3 GeV ||40|]; SLAG Mark I data from 3 to 6 GeV [y]; and Grystal Bah data above 5 



GeV [42|. For useful compilations and reviews see Ref. [43|. We used simple fits to the data 
in some regions, particularly when the data had a lot of structure and/or had large statistical 
errors. This was more convenient for the numerical integration routine and made it easier for 
us to examine the effect of the experimental uncertainties on the smeared result. Fig. 4 shows 
our data compilation, up to 6 GeV, excluding narrow resonances. In fact, we used data going 
well beyond h threshold, but they have no real effect on the results we present. 
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The sharp resonances uj, <j), J/tp, and '(/'(3770) were not included in the data compilation 
so that their contribution to Rpqw could be put in analytically. They have a relativistic Breit- 



Wigner form |44]: 

9 M^r^ 

where M, F, Bu, and are, respectively, the mass, width, leptonic branching fraction, and 
hadronic branching fraction of the resonance. The parameters for the resonances were taken 
from the 1992 Review of Particle Properties |^6|. For Bu we used the weighted average of the 
ee and /i/x branching ratios. 

The contribution of such a Breit-Wigner resonance to the smearing integral ( p. 45 ) can be 
evaluated analytically using partial fractions. [The resulting expression is too cumbersome to 
quote, but we may note that the narrow width approximation: 

1 vr 



(s - M2)2 + M2F2 ^ MF 
which gives a contribution to Rpqw of 



5{s - M2), (3.48) 



- 9BiiBhAMT 

is a pretty good approximation.] In Fig. 5 we show, for two different A values, the contributions 
of the various resonances to the experimental Rpqw ■ The p's contribution is shown by a dotted 
line. However, since the p is rather wide and asymmetric, it was actually treated, not in this 
manner, but by numerical integration, using the data points from Ref. |3^ as part of the data 
compilation (Fig. 4). 



3.4 Results and uncertainty estimates 

The results obtained by applying PQW smearing to both theory and experiment are shown in 
Fig. 6. For the smaller A (1 GeV^) there is good agreement between theory and experiment 
below 1 GeV, but in the charm-threshold region there is clearly insufficient smearing for the 
comparison to be meaningful. Increasing A to 3 GeV^ smooths out the experimental curve 
almost completely. The agreement between theory and experiment is excellent below 2 GeV. In 
the charm region the agreement is less good, but this can be attributed mainly to the sizable 
systematic normalization uncertainty (10 - 20%) in the data in this region, which produces an 
uncertainty of about ±0.4 in the experimental Rpqw at around Q = 4 GeV. For comparison. 
Fig. 6(b) also includes the naive parton-model prediction. One can see from this that the QCD 
correction term TZ provides about a 20% increase which is vital to the good agreement with the 
data. 
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Using Q^, rather than Q as the variable, we can continue RpqwiQ"^) into the negative 
region (Cf. Ref. [^Bj). As shown in Fig. 7, for A = 1 GeV^, the good agreement persists. 

To quantify the good agreement at low energies, we discuss how various uncertainties would 
affect RpQw at Q = 0. First we discuss the experimental uncertainties. There is about a 
5% uncertainty in the p, uj and (j) contributions, due to the uncertainty in their total and 
leptonic widths. For A = 1 (3) GeV^ this gives an error in RpQw{0) of about ±0.04 (±0.02). 
Uncertainties in the Tp resonance parameters affect RpQw{0) by ±0.01 or less. We considered 
the effect of a 15% normalization change in the continuum data in the 1.5 - 3 GeV region: The 
effect on Rpqw{0) was about ±0.03 (±0.06) for A = 1 (3) GeV^. We also allowed for a 15% 
normalization change in the 3-5 GeV region. The effect on ^pqvk(O) was about ±0.02 (±0.065) 
for A = 1 (3) GeV^. Combining these four distinct sources of error in quadrature, we estimate 
an overall uncertainty in the experimentally determined RpQ^iO) of ±0.06 for A = 1 GeV^ and 
±0.08 for A = 3 GeV^. 

On the theoretical side, errors arise from two sources: (i) uncertainty in the input parameters 
(quark masses and Aj^), and (ii) truncation of the perturbation series. We varied each quark 
mass by its quoted error |^^. Varying the u and d masses had negligible effect. RpQw{0) changed 
by ±0.004 (±0.001) on varying the s mass, and by ±0.005 (±0.013) on varying the c mass for 
A = 1 (3) GeV^. Changing Aj^ by 50 MeV to 280 MeV increased Rpqw{0) by 0.019 (0.014) 
for A = 1 (3) GeV^. The series-truncation error can reasonably be estimated from the last 
term in the optimized series, as we argued earlier. At 1 GeV this suggests that TZ is accurate 
to about 10%, and is considerably more accurate at larger energies. This is corroborated by 
the good agreement between second- and third-order results. The theoretical uncertainty in TZ 
above 1 GeV contributes an error in Rpqw{0) of less than ±0.006 (±0.009) for A = 1 (3) GeV^. 
Below 1 GeV the prediction for TZ is much more uncertain. However, as discussed in Subsect. 
2.6, we think that even at Q = the result is reliable to within a factor of 2. Conservatively, 
we considered the effect of increasing the predicted 7^ by a factor of 2 over the whole range, 
< Q < 1 GeV. This affects the low-energy Rpg^,- by 0.033 (0.011) for A = 1 (3) GeV^. If we 
linearly add all the above-mentioned uncertainties we get a total uncertainty of ±0.07 for A = 1 
GeV^ and ±0.05 for A = 3 GeV^. Thus, the theoretical uncertainties are comparable to the 
experimental uncertainties. 

3.5 Significance of the results 

We can now discuss the significance of the agreement between theory and experiment. We first 
ask: How restrictive is the data? To quantify the discussion we define a 'straw-man' model for 
TZ in which TZ is the same as the OPT result down to 2 GeV, but then follows the one-loop. 
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3-flavour form, (12/27)(l/ln Q^/Ag), with Aq ~ 0.2 GcV, until it reaches a value H, at which 
it remains frozen down to Q = 0. If the "freeze-out" value, H, is about 0.3, then this 'H model' 
is essentially equivalent to the OPT result. If H is much larger then this model gives a result 
for RpQw{0) that is too large by more than the uncertainties just estimated. We find that H^s 
above 2 are disfavoured by the data. (As an illustration Fig. 8 shows the result with H = 4.6, 
which is clearly ruled out.) At the other extreme, the data disfavour an H less than 0.09. Thus, 
although a wide range of H values can be tolerated, the data do imply that the couplant cannot 
grow very large in the infrared region; nor can it remain too small. 

Next we ask: How predictive is the theory? Because of the need for smearing, the theory 
tells us almost nothing about the shape or structure of the e+e~ data in the region below 1 
GeV. However, it does tell us something about the average magnitude of the cross section. 
The low-energy data is, in fact, dominated by the p peak. After smearing with A = 1 GeV^, 
this contributes about 0.7 to Rpqw below 1 GeV. Thus a 10% change in the area under the p 
peak would change Rpqw by the ±0.07 estimated uncertainty in the theoretical prediction. We 
conclude that perturbative QCD can tell us, at least crudely, the size of the p resonance. 

3.6 The smeaired derivative 

As an extension of PQW's ideas we also considered a quantity: 



This provides a somewhat different test, though obviously not an independent test, of the rela- 
tionship between theory and experiment. Its calculation requires only straightforward modifi- 
cations to the procedures used to calculate Rpqw 

In Fig. 8 we compare the smeared derivatives from theory and experiment for A = 2 and 4 



gives the first peak just below 3 GeV. However, there is clearly insufficient smearing in the 
charm region. Increasing A to 4 GeV^ greatly smooths out both curves and gives quite good 
agreement. 

4 Phenomenological Virtues of a Frozen Couplant 

The idea that the strong coupling constant, as{Q'^), "freezes" at low energies has long been a 
popular and successful phenomenological hypothesis. We first note that a freezing of as{Q'^) is a 




(3.50) 



which represents a "smeared derivative" , in the sense that 



limD(Q,A) = di?e+e-/dg'- 



(3.51) 



Gev^. For A = 2 GeV^ there is good agreement at low energies, and the theory qualitatively 
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natural consequence of a picture where the gluon aquires an effective, dynamical mass rUg 
Naively, this would modify the leading-order, 3-flavour couplant to: 

TT 27 ln[(Q2 + 4m2)/A2] ' ^ ^ 

a form that has been used in many phenomenological papers. For mg a little larger than Aq 
this gives a zero-Q value comparable to ours. Note, though, that the variation with Q at low 
energies is somewhat different from ours in Fig. 2. 

Another commonly used form is the "hard-freeze" form in which: 



a,(g2) J (12/27)(l/lnQVA2), for > 

(4.53) 



TT 



constant = H, for < Qq, 



with H = (12/27)(l/lnQg/Ag). This is the "i? model" that we mentioned in Subsect. 3.5. For 
H 0.26 {i.e., Qq/Kq « 2.3) it is a close approximation to our as{Q'^)/T^ shown in Fig. 2. 

We now briefly survey some of the phenomenological literature in order to make two points: 

(i) a frozen Us provides a way to understand many important facts in hadronic physics, and 

(ii) the values extracted phenomenologically are very much in accord with our low-Q value 
as/vr = 0.26. [Note that we quote as/vr rather than values.] 

(a) Total hadron-hadron cross sections, although slowly rising at very high energies, are 
remarkably constant over a wide energy range, and their relative sizes correlate with their quark 
content in a very suggestive way. A simple and succesful description is provided by the 2-gluon- 



exchange model |47], based on the Low-Nussinov model of the Pomeron |^^. This model requires 



a finite couplant at low momentum transfer, and Ref. |47| found a value Os/vr ~ 0.17. A recent 
version of this model, framed in terms of a dynamical gluon mass {mg = 0.37 GeV, for Aq = 0.3 
GeV), is given in Ref. |49]. Another recent version of this model uses the 'H-model' form of 
0(s{Q'^)/'^- In order to fit the absolute magnitude of the vr-nucleon cross section, Qq needs to be 
about 0.44 GeV if Aq = 0.2 GeV. This corresponds to H = 0.28. The same frozen couplant 
has been used successfully in subsequent work on deriving nucleon structure functions from the 
constituent-quark model | |5l[ |. 

(b) Hadron spectroscopy also points to a low-energy couplant of around 0.2 ~ 0.25 
Godfrey and Isgur ||5^ provide a unified description of light- and heavy-meson properties in a 
"relativized" potential model with a universal one-gluon-exchange-plus-linear-confinement po- 
tential. For the model to work for light mesons it is crucial to incorporate relativistic effects, 
and to employ a form of the couplant that freezes at low energies. Their fits yield a form of 
o.s{Q'^)/'^ that freezes to about 0.19, and has a shape similar to ours. In a fully relativistic 



treatment Zhang and Koniuk [54| can naturally explain why the vr is so much lighter than the 
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p. The tt/ p mass ratio is a steeply falling function of the strong couplant, and the experimental 
value occurs at Os/vr = 0.265 

(c) Hadron form factors at low energies can be successfully treated assuming a frozen cou- 
plant, as shown in Ref. |55], which used the form ( |4.52| ) with nig ~ 0.1 to 0.5 GeV, for Aq « 0.1 
GeV. 

(d) Chiral soliton models of the nucleon can fit a wide variety of nucleon properties if one 



includes one-gluon exchange corrections with an a^/vr of about 0.2 |53|. Ref. finds that 
the experimental deviation from the Gottfried sum rule, the A-nucleon mass difference, the first 
moment of the polarized proton structure function, and the neutron-proton mass difference all 
require a common as/vr value. (However, the actual value found, 0.2, could be re-scaled by 
making a different choice for another parameter in the model 15^.) Other nucleon properties 



are consistent with an as/vr of this size [57]. 

(e) The pT spectrum in W, Z production in pp or pp collisions can be successfully predicted 
by QCD right down to px = if multiple gluon radiation effects are appropriately re-summed 
|l58| . However, it is essential in the low-px region to invoke a freezing of as(p^). The form ( [4.52| ) 
has been used, with 4m^/AQ denoted by 'a'. Unfortunately, the results are very insensitive to 



the parameter a; anything in the range 3 - 100 gives an acceptable fit to current data [59 



This corresponds to a range 0.1 - 0.4 for the zero-Q couplant. Perhaps, future data will make 
it possible to narrow this range. 

(f) Jet properties can be quite successfully described by the "modified leading log approxi- 
mation" |6^ , but to obtain predictions at small momenta it is necessary to invoke a freezing 
of the couplant. Fits to data on heavy-quark-initiated jets give zero-Q values of Os/vr around 
0.22 [^]. This value depends somewhat on the form of OsiQ"^) assumed in the fit, but it was 
found empirically that the result for the integral: 

/ d/fc ^^^^ « 0.2 GeV (4.54) 
Jo vr 

was fit invariant [^]. Integrating our a^/vr in Fig. 2 leads to precisely 0.2 GeV. 

(g) Hadron-hadron scattering at very high energies where the cross sections rise asymptoti- 



cally, but must satisfy unitarity, seems to call for the 'critical Pomeron' picture |63], at least as 
a first approximation. It seems that one could only hope to derive such a picture from QCD if 
there is an infrared fixed point [^4|. In fact. White has argued for additional quarks, or colour- 
sextet quarks, in order to have Nf effectively equal to 16 |Q. (See the discussion in Subsect. 
2.6 above.) However, our results imply that the infrared fixed point persists down to low Nf. 
This may mean that one can have all the virtues of White's picture without the need for more 
quarks. 
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5 Summary and Conclusions 



We have applied OPT to the third-order QCD calculation of i2g+g-. At energies above about 
1 GeV there is every sign that the approximation is healthy: the perturbation series in the 
"optimized" scheme is well behaved, and there is good agreement between second- and third- 
order results. This was not true of the situation created by the old, incorrect R^+e- calculation 



, 31, ^] (see Table 2). The contrast between the 'old' and 'new' results emphasizes the point 



that the third-order i2e+e- calculation provides a very real, empirical test of "optimization" 



ideas. At the time, the statements [32] that the 'old' third-order results [|l| tended to cast doubt 
on the usefulness of "optimization" were perfectly fair comment. Because of this history we take 
especial satisfaction in the transformed situation produced by the new result 

Furthermore, contrary to the old situation, the optimized couplant now shows a marked 
decrease from second to third order. This is in accord with the "induced convergence" conjecture 
that "optimization" naturally cures the divergent-series problem ||5|, [l^ . 

The third-order OPT results remain finite down to Q = 0, with the optimized couplant, 
Os/tt, "freezing" to a value 0.26 below 300 MeV. No ad hoc assumption was used to obtain this 
result: it is the direct consequence of using the calculated R^+e- and j3 series coefficients as 
inputs to the "optimization" procedure specified in Ref. Q. 

It must be admitted that, at very low energies, the prediction for IZ (the QCD correction 
term in R^+^-) has a large uncertainty. Since third order is the lowest order at which it is 
even possible to get finite infrared results, one should not be surprised if the approximation is 
somewhat crude. Nevertheless, as we discussed in Subsect. 2.6, the qualitative conclusion that 
TZ remains small (say, 0.3 it 0.3 at Q = 0) is inescapable in the context of OPT. 

The OPT prediction is supported by the data. As we showed in Sect. 3, the PQW-smeared 
Re+e- data is consistent with a perturbative QCD description, provided that the couplant freezes 
to a modest value at low energies. 

The hypothesis that the couplant freezes at low energies has been used very successfully in 
a wide variety of phenomenological work, where the low-energy couplant is treated as a free 
parameter to be fitted to experiment. The values that emerge are quite comparable to ours. 
There are some other theoretical indications of a freezing of the couplant |4^, but our 
evidence is remarkable in that it comes solely from perturbation theory and RG invariance. The 
predicted value, Os/tt = 0.26, for the frozen couplant is a purely theoretical number. It does not 
depend on knowing the value of Aj^, but only on knowing the number of light quarks. 
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Appendix A: Numerical solution of the optimization equations 



After expressing ri and r2 in terms of a and C2, using the pi, p2 definitions and ( p.42 ), we have to 



simultaneously solve the optimization equations, ( |2.27 , 2.28| ). These define two curves in the a 



C2 plane whose intersection point we seek. In what we call the 'spiralling' method, ( 2.2g| ) is first 
solved for a; then, with this a, ( |2.27| ) is solved for C2; then, with this C2, ( |2.28| ) is solved for a; 
and so on. Which equation is solved for which variable is crucial; the other choice would 'spiral 
out' from the desired solution. (The standard "secant method" |6(:] was generally sufficient for 
solving the individual optimization equations.) 

A convenient starting point for this iterative procedure was provided by an approximate solu- 
tion to the optimized equations due to Pennington, Wrigley, and Mignaco and Roditi (PWMR) 
|l67|| . This approximation expands the optimization equations ( |2.27| , ^.281) as a series in a and 
keeps only the lowest non-trivial term. Noting that / = a(l — ca + ...), it is easy to check that 
this gives: 

fi -3C2. (5.55) 

One can improve this approximation by writing each fj as a series in a and successively equating 
coefficients of different orders in a to zero. To next order this gives: 

fi ic2a, f2 ^ -ic2 + icc2a. (5.56) 

[One may note that in TZ^^^ there is a near cancellation between the second and third order 
terms, fia and f2a^. Thus, IZ^^^ turns out to be closely equal to a^^\] 

The 'spiralling' method worked well for Q > 0.3 GeV starting from the PWMR solution. 
However, at lower energies the PWMR approximation breaks down, and does not provide a 
satisfactory initial guess. In fact, at the infrared fixed point one has instead, from ( p. 35 ): 



n = ic2a*, r2 = -§C2. (5.57) 

We therefore proceeded to low Q in successive stages, utilizing the solution at the previous Q 
as the initial guess for the next lower Q. We also encountered a 'creep' problem: At low Q the 
two curves representing the "optimization" equations become almost parallel (each being almost 
parallel to the infrared boundary line 1 + ca + 020^ = 0) and they cross at a very small angle. 
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Thus, instead of 'spiralling in' to the solution, one creeps towards it stepwise. The convergence is 
very slow and the danger is that the solution can appear to have converged within the specified 
tolerance, when in fact it still has a considerable way to go. To avoid this pitfall we would repeat 
the procedure from a different starting point, so as to creep towards the solution from the other 
side. In this way we could bracket the true solution, and hence ensure reliable accuracy. 

We also tried the 'intersection' method as an alternative. Taking an initial guess for C2, one 
solves for a in each of the two optimized equations. For each a one then solves the other equation 
for C2- This gives a pair of points on each of the two curves. The straight lines that join up each 
pair should approximate the curves themselves, and hence their intersection should approximate 
the desired solution. The procedure can then be iterated. This method also worked well for 
Q > 0.3 GeV starting from the PWMR solution. At lower energies, where the two curves become 
nearly parallel, this method did not suffer from the 'creep' problem, but it had the opposite vice: 
it tended to make such a large extrapolation in each iteration that it would become unstable 
and erratic. 

Appendix B: Flavour thresholds 

Since TZ has been calculated only with massless quarks, we are really approximating "full QCD" 
with a set of effective theories, each with a different number of massless quarks. The 
parameters of these theories need to be appropriately matched, so that they correspond to 
a single, underlying "full QCD" theory. The point is well explained by Marciano |^8|, who 
provides explicit formulas for matching Aj^ across thresholds. Unfortunately his analysis uses 
a truncated expansion of a{ij,) in powers of 1/ ln(^/A), which would not be a valid approximation 
at low energies; in particular at s-quark threshold. 

Our procedure was simply to require the optimized TZ^^^ to be continuous at a threshold. This 
was done numerically by running our optimization program at the threshold energy (Q = 2mq) 
with both values of Nf and adjusting one of the A^ parameters until the two 'RS^'^ results agreed. 
Starting with = 230 MeV for 4 fiavours, we found = 281 MeV, and A^^ = 255 MeV. [In 

MS ' MS ' MS ^ 

terms of A the corresponding values are: 257, 308, and 277 MeV for 4, 3, 2 flavours, respectively.] 
Essentially the same results were obtained if we required instead that a be continuous. We 
checked that this procedure agreed very closely with Marciano's formulas at both c- and 6-quark 
thresholds. 

It is noteworthy that we find Ai^i. to be smaller than Ai^^, contrary to the pattern at the 
higher thresholds. Our final results are very insensitive to the A^ value, however, because at 
energies below s threshold the TZ results are essentially governed by the infrared fixed point. 
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AT 


P2 







-8.410 


0.313 


1 


-9.997 


0.280 


2 


-10.911 


0.263 


3 


-12.207 


0.244 


4 


-13.910 


0.224 


5 


-15.492 


0.208 


6 


-17.665 


0.191 



Table 1: p2 invariants and fixed-point couplants for Nf = to 6. 





5, Q = 34 GeV 


Order 


a 




r2 


7^ 


change 






2nd 


0.0415 


-0.599 




0.0404 






= 100 MeV 


3rd 


0.0394 


-0.301 


7.64 


0.0394 


-2.4% 




= 87 MeV) 


old 


0.0452 


+1.363 


-29.48 


0.0453 


12% 






2nd 


0.0569 


-0.588 




0.0550 




^M5 ' 


= 500 MeV 


3rd 


0.0526 


-0.405 


7.71 


0.0526 


-4.4% 




= 436 MeV) 


old 


0.0690 


+1.988 


-27.59 


0.0694 


26% 



Table 2: Comparison of second- and third-order optimized results: 'old' refers to third order 
with the old, incorrect result. 
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Q (GeV) 


Nf 


a 


ri 


r2 


7^ 


error 


3.0 


4 


0.076 


-0.53 


6.9 


0.076 


4% 


1.0 


3 


0.126 


-0.79 


6.3 


0.126 


10% 


0.4 


3 


0.221 


-1.77 


8.8 


0.229 


43% 





2 


0.263 


-2.89 


14.6 


0.330 


100% 



Table 3: Illustrative third-order optimized results at low energies. Ajjg{4 flavours) = 230 MeV. 
The estimated fractional error is |f2a^|. 
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Figure Captions 



Fig. 1. The optimized solutions in the a, C2 plane for 2 and 3 quark flavours, in the low-energy 
region. The open squares represent the fixed-point solution, Eqs. (2.37, 2.36), which lies on 
the infrared boundary 1 + ca + C2a^ = 0. The boundary is shown by the solid line (Nj = 2) 
or the dashed line {Nj = 3) at the right. The dotted vertical lines are to indicate ss threshold 
at Q = 0.40 GeV where Nf changes from 3 to 2. (The A parameters are matched so that R is 
continuous (see Appendix B.), but there are then slight discontinuities in a and C2.) The dotted 
line shows the solution for a 3 flavor world down to Q = 0, while the dashed line shows a 2 
flavour world extending up towards 1 GeV. The points shown are spaced at 0.05 GeV intervals 
from Q = 0.40 GeV. 

Fig. 2. The optimized third-order results for a = as /it and 1Z^^\ Also shown is the second- 
order result, iZ^'^\ Quark thresholds are indicated by the vertical lines. 

Fig. 3. The perturbative QCD prediction for R^+e- from third-order OPT (solid line). The 
inset shows the region around u and d quark thresholds. The dashed line is the parton- model 
prediction. 

Fig. 4. Compilation of experimental Rf,+e- data (excluding narrow resonances). A few 
representative statistical error bars are shown. The solid line represents an 'eyeball fit'. 

Fig. 5. The contributions of narrow resonances to Rpqw for two values of the smearing 
parameter A. 

Fig. 6. Comparison of "smeared" theoretical and experimental results. The parton-model 
result is shown by the dotted line in (b). 

Fig. 7. Comparison of "smeared" results extended to spacelike Q^. The dotted line shows a 
'straw-man' model in which the couplant becomes large at low energies (see subsect. 3.5). 

Fig. 8. Comparison of theoretical and experimental results for the "smeared derivative" (Eq. 
( 3.50| )) for two values of A. 
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Fig. 1. The optimized solutions in the a, C2 plane for 2 and 3 quark flavours, 
in the low-energy region. The open squares represent the flxed-point solution, Eqs. 
(2.37,2.36), which lies on the infrared boundary 1 + ca + c^a^ = 0. The boundary 
is shown by the solid line (iVj = 2) or the dashed line (iVj = 3) at the right. The 
dotted vertical lines are to indicate ss threshold at Q = 0.40 GeV where Nj changes 
from 3 to 2. (The A parameters are matched so that R is continuous (see Appendix 
B.), but there are then slight discontinuities in d and C2.) The dotted line shows 
the solution for a 3 flavor world down to Q = 0, while the dashed line shows a 2 
flavour world extending up towards 1 GeV. The points shown are spaced at 0.05 
GeV intervals from Q = 0.40 GeV. 




Fig. 3. The perturbative QCD prediction for Rg+e- from third-order OPT (sohd 
hne). The inset shows the region around u and d quark thresholds. The dashed 
hne is the parton-model prediction. 




Fig. 4. Compilation of experimental Rg+e- data (excluding narrow resonances). 
A few representative statistical error bars are shown. The solid line represents an 
'eyeball fit'. 




Fig. 7. Comparison of "smeared" results extended to spacelike . The dotted 
line shows a 'straw-man' model in which the couplant becomes large at low energies 
(see subsect. 3.5). 




Fig. 2. The optimized third-order resuhs for a = as/n and T^^"^^ Also shown is the 
second-order resuh, TZ^^K Quark thresholds are indicated by the vertical lines. 




Fig. 6. Comparison of "smeared" theoretical and experimental results. The 
parton-model result is shown by the dotted line in (b). 
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Fig. 5. The contributions of narrow resonances to Rpqw for two values of the 
smearing parameter A. 




Fig. 8. Comparison of theoretical and experimental results for the "smeared 
derivative" (Eq. 3.50) for two values of A. 



